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1 Introduction 


The predictions of the standard model (SM) JT}j7| of fundamental interactions have been con¬ 
firmed by a large number of experimental measurements. The observation of a new boson 
with a mass of 125GeV and properties compatible with those of the SM Higgs boson HSl-fTOl. 
confirms the mechanism of the electro weak symmetry breaking (EWSB). Despite the success of 
this theory in describing the phenomenology of particle physics at present collider energies, the 
mass of the Higgs boson in the SM is not protected against quadratically divergent quantum- 
loop corrections at high energy Supersymmetry (SUSY) IflTl 12] is one example of alternative 
models that address this problem. In SUSY, such divergences are cancelled by introducing a 
symmetry between fundamental bosons and fermions. 

The minimal supersymmetric extension of the standard model (MSSM) Ifl3l fl4| predicts the 
existence of two Higgs doublet fields. One doublet couples to up-type and one to down-type 
fermions. After EWSB, five physical Higgs bosons remain: a CP-odd neutral scalar A, two 
charged scalars H^, and two CP-even neutral scalar particles h and H. The neutral bosons 
h. A, and H, will be generically referred to as c[> collectively in this paper, unless differently 
specified. 


At lowest order in perturbation theory, the Higgs sector in the MSSM can be described in terms 
of two free parameters: m,\, the mass of the neutral pseudoscalar A, and tan /3, the ratio of the 
vacuum expectation values of the two Higgs doublets. The masses of the other four Higgs 
bosons can be expressed in terms of these two parameters and other measured quantities, such 
as the masses «w and niz of the W and Z bosons, respectively. In particular, the masses of the 
neutral MSSM scalar Higgs bosons H and h are given Itl3l by 


«H,h = 


\ {m Z A + m Z ± [( m A + m z) 2 - 4»d ? ”z c °s 2 1/2 } 


i 1/2 


( 1 ) 


The A and H bosons are degenerate in mass above 140 GeV and for small cos (5 (large tan /3) 
values. This expression also provides an upper bound on the mass of the light scalar Higgs 
boson, corresponding to nih < mz|cos2j3|. The value can become as large as « 135 GeV 
once radiative corrections are taken into account Itl5l . 


The main production mechanisms for the three neutral cp bosons at the LHC are the associated 
production with bb quarks (AP), given at the leading order by the Feynman diagram shown in 
Fig. [l](left), and the gluon fusion (GF) process, shown in Fig. [T|(right) IflhHlBl . The GF process 
with virtual t or b quarks in the loop is dominant at small and moderate values of tan ft. At 
large tan /I the coupling of tp to down-type quarks is enhanced relative to the SM Ifl9l and the 
AP process becomes dominant. Similarly, the coupling of the (p boson to charged leptons is also 
enhanced at large tan /!. 

This paper reports on a search for the MSSM neutral Higgs bosons produced either by the AP 
or GF mechanisms, where the Higgs bosons decay via cp —> ji l p ■ The analysis is sensitive 
to all the three bosons, h, H, and A in the mass range between 115 and 300 GeV. The search 
is performed by the CMS collaboration using data recorded in pp collisions at the LHC, cor¬ 
responding to an integrated luminosity of 5.1 fb 1 at y/s = 7TeV and 19.3 fb 1 at ^/s = 8TeV. 
The common experimental signature of the two processes is a pair of oppositely charged muons 
with high transverse momentum ( pj ) and a small imbalance of pj in the event. The AP process 
is characterized by the presence of additional jets originating from b quarks (b jets), whereas the 
events with only jets from light quarks or gluons are sensitive to the GF production mechanism. 
The presence of a signal would be characterized by an excess of events over the background in 
the dimuon invariant mass corresponding to the (p mass value. 
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2 The CMS detector and event reconstruction 




h, H, A 


Figure 1: Leading-order diagrams for the main production processes of MSSM Higgs bosons at 
the LHC (left) in association with bb production and (right) through gluon fusion. 

Although the product of the cross section and the branching fraction for the fi 1 channel is 
a factor 10 3 smaller than for the corresponding final state, the muon pair can be fully 

reconstructed, and the invariant mass precisely measured by exploiting the excellent muon 
momentum resolution of the CMS detector. Searches for the MSSM Higgs bosons have been 
performed at LHC by the LHCb experiment in the t - t final state at large pseudorapidity 
values m, the ATLAS experiment in the }t 1 }i and t + t channels |TT[[ 22 ] , and by the CMS 
experiment in the r 1 r l23l and bb EH [251 final states. Limits on the existence of MSSM Higgs 
bosons were also determined at Tevatron E6ti29l and at LEP l30l . 

Traditionally, searches for MSSM Higgs bosons are presented in the context of benchmark 
scenarios that describe the mass relation among the three neutral MSSM Higgs bosons, their 
widths, and cross sections. Each scenario assigns well defined values to the relevant param¬ 
eters of the MSSM, except m,\ and tan /3, which are left free to vary. The m™ ax benchmark 
scenario EH [3T| provides m h values as large as 135 GeV, and the weakest bounds on tan /l for 
fixed values of the top quark mass. For this reason, it has been used in most of the previously 
quoted analyses to present the results from MSSM Higgs boson searches. However, within 
the MSSM the newly discovered state with a mass of 125 GeV can be interpreted as the light 
CP-even Higgs boson, h l32l . In this case, a large part of the m A -tan fi parameter space is ex¬ 
cluded within the m™ ax scenario, and new benchmarks were therefore proposed in which the 
MSSM parameters are adjusted to have m h in the interval 122 to 128 GeV, but with a wider 
range of tan fi and rn,\ values EH [3Tj [32l. To do this, the m™ ax scenario was reformulated in 
two versions, m™ od+ and »zj(' od corresponding to different values of the top squark mixing 
parameter. Other recently proposed scenarios 1311 are the light top squark (light stop) model, 
which results in a modified GF rate, and the light tau slepton (light stau) model, which yields 
a modified h —» 77 branching fraction. Such models are expected mainly to affect the Higgs 
boson production cross section and not the kinematic properties of the events. A list of the pa¬ 
rameters of the various scenarios can be found in Ref. E3l . The results presented in this paper 
are obtained in the framework of the MSSM m™ d+ scenario. Comparisons are also made with 
other benchmarks. 


2 The CMS detector and event reconstruction 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintilla¬ 
tor hadron calorimeter (HCAL), each comprised of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Forward calorimetry extends the coverage provided by the barrel and endcap detec- 
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tors up to pseudorapidity |)/1 < 5. A detailed description of the CMS detector, together with 
a definition of the coordinate system and kinematic variables, can be found in Ref. Il33l . The 
CMS offline event reconstruction creates a global event description using the particle flow (PF) 
technique |34| . The PF event reconstruction attempts to reconstruct and identify each particle 
with an optimized combination of all subdetector information. The missing pj vector is de¬ 
fined as the projection on the plane perpendicular to the beams of the negative vector sum of 
the momenta of all reconstructed particles in an event. Its magnitude is referred to as E™ ISS . 

An average of 9 and 21 pp collisions take place in any LHC bunch crossing, respectively at 7 
and 8 TeV, because of the large luminosity of the machine and the size of the total inelastic cross 
section. These overlapping events (pileup) are characterized by small-pj tracks, compared to 
the particles produced in a <p —> }i + p~ event, and their presence can degrade the detector 
capability to reconstruct the objects relevant for this analysis. The primary vertex is chosen 
from all reconstructed interaction vertices as the one with the largest sum in the squares of 
the pj of the associated tracks. The charged tracks originating from another vertex are then 
removed. 

Offline jet reconstruction is performed using the anti-fcj clustering algorithm [35], 36] with a dis¬ 
tance parameter of 0.5. The jet momentum is defined by the vectorial sum of all the PF particles 
momenta in the jet, and found in simulation to be within 5% to 10% of the true hadron-level 
momentum, with some pj and ;/ dependence. Extra energy coming from pileup interactions 
affects the momentum measurement. Corrections to the measured jet energy are therefore ap¬ 
plied. They are derived from event simulation, and confirmed with in-situ measurements using 
energy balance in dijet and Z/photon+jet events |37|. 

Muons are measured in the pseudorapidity range \rj\ < 2.4, using detection planes based on 
three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers. Matching 
muons to tracks measured in the silicon tracker provides relative pj resolutions for muons 
with 20 < pj < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcaps. The pj 
resolution in the barrel is better than 10% for muons with pj up to 1 TeV l38l . 


3 Simulated samples 

Simulated samples are used to model the signal and to determine the efficiency of the signal 
selection. Background samples are also simulated to optimize the selection criteria. The nor¬ 
malization and distribution of the background events are measured from data. 

The signal samples are generated using the Monte Carlo (MC) event generator pythia 6.424 Il39l 
for a wide range of m a and tan /3 values, as listed in Table [l] for the AP and the GF produc¬ 
tion mechanisms. The <p production cross sections and their corresponding uncertainties are 
provided by the LHC Higgs Cross Section Working Group I16W181 . The cross sections for the 
GF process in the m“ ax scenario are obtained using the HIGLU program f40l [41], based on 
next-to-leading order (NLO) quantum chromodynamics (QCD) calculations. The SUSHI pro¬ 
gram C42] is used for the other benchmarks. For the AP process, the four-flavor NLO QCD 
calculation Il43l [44] and the five-flavor next-to-next-to-leading order (NNLO) QCD calculation 
are implemented in BBH@NNLO |45] and combined using the Santander matching scheme 06]. 
The Higgs Yukawa couplings computed with the FEYNHIGGS program 14 7\ are used in the cal¬ 
culations. The decay branching fractions to muons in the different benchmark scenarios are 
obtained with FEYNHIGGS and HDECAY [|48| . Further details on signal generation can be found 
in Refs. OME]. 
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4 Event selection 


The values of predicted by FEYNHIGGS differ typically by a few GeV from those computed 
with PYTHIA. The invariant mass spectrum of the h boson is therefore shifted to match the 
FEYNHIGGS prediction. The small difference between PYTHIA and FEYNHIGGS in assessing the 
width of the h boson is of the order of 100 MeV, and therefore neglected, since the experimental 
mass resolution is at least one order of magnitude larger. The PYTHIA parameters used to sim¬ 
ulate the signal are those for the m™ ax scenario. Since for a given set of m,\ and tan /3 values, the 
kinematic properties of the final state are the same for all the scenarios, the simulated samples 
based on the m™ ax benchmark are also used to check the validity of the other models. Further 
details on this procedure and the related systematic uncertainties are discussed in Section [7] 

Table 1: The Ma and tan/I values used to generate signal samples. 


m A (GeV) 

hza step (GeV) 

tan /3 

tan j6 step 

115-200 

5 

5-55 

5 

200-300 

25 

5-55 

5 

300-500 

50 

5-55 

5 


The main source of background for the (p production and decay to p 1 }i is Drell-Yan muon- 
pair production, qq —I Z/y* —I /by . Another background is from oppositely charged muon 
pairs produced in decays of top quarks in tt production. These events are simulated using the 
MadGraph 5.1 |49l generator. Other background processes such as W^W - , W ± Z, and ZZ are 
generated with PYTHIA. The MC samples also include simulated pileup events to reproduce 
the overlapping pp interactions present in the data. All generated events are processed through 
a detailed simulation of the CMS detector based on Geant4 H50l and are reconstructed with 
the same algorithms used for data. 


4 Event selection 

The experimental signature of the MSSM Higgs bosons decay considered in this analysis is a 
pair of oppositely charged muons with high pj. The invariant mass of the pair corresponds to 
the mass of the (p boson within the experimental resolution. Moreover, the process is character¬ 
ized by a small Ey" ss in the event. If the (p boson is produced in association with a bb pair, the 
presence of at least one b quark jet is expected. 

The details of the event selection are listed below, and summarized in Table[2] The events are se¬ 
lected using a single-muon trigger, which requires at least one isolated muon with pj > 24 GeV 
in the pseudorapidity range \rj \ < 2.1. The distance of the primary vertex along the z axis from 
the nominal centre of the detector must be |zpy| < 24 cm. Muon candidates are reconstructed 
and identified using both the inner tracker and the muon detector information. The selected 
events must have at least two oppositely-charged muon candidates, each with pj > 25 GeV. 
In events with more than two muon candidates, the two with opposite charges and the high¬ 
est pj are retained. The rj of both muon candidates is chosen to match the trigger acceptance. 
Each muon track must have at least one hit in the pixel detector, more than five or eight lay¬ 
ers with hits in the tracker, respectively, for the 8 and 7 TeV data and a directional matching 
to hits in at least two different muon detector planes. In addition the global fit to the hits of 
the muon candidate must include at least one hit in the muon detector. The y 2 / dof of the 
global fit of the muon track must be smaller than 10. These requirements ensure a good mea¬ 
surement of the momentum, and significantly reduce the amount of hadronic punch-through 
background H38l . To reject cosmic ray muons, the transverse and longitudinal impact param¬ 
eters of each muon track must satisfy the requirements \d xy \ < 0.02cm and \d z \ < 0.1cm, 
respectively. Both parameters are defined relative to the primary vertex. To ensure that the trig- 
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Figure 2: The Ey ISS distribution for events with a reconstructed dimuon invariant mass 
rn^ i„ > 60GeV in data and in simulated events at y/s = 7 (left) and \fs = 8TeV (right). 

The expected contribution is also shown for a signal at m f \ = 150 GeV and tan /3 = 30. 



Figure 3: The distribution of the b tagging discriminant, b^ isc , for events that satisfy the selection 
Ey ISS < 35 GeV in data collected at y/s = 7 (left) and yfs = 8TeV (right). For each event, the 
largest value of &disc is selected. 

ger muon candidate is well-matched to the reconstructed muon track, at least one of the two 
muon tracks is required to match the direction of the trigger candidate within a cone A R = 0.2, 
where A R = V(A?/) 2 + (A q>) 2 is the distance between the muon track and the trigger candi¬ 
date direction in the t]-cp plane, with cp being the azimuthal angle measured in radians. Both 
reconstructed muon candidates must fulfill isolation criteria. A muon isolation variable is con¬ 
structed using the scalar sum of the px of all PF particles, except the muon, reconstructed within 
a cone A R = 0.4 around the muon direction. A correction is applied to account for the possible 
contamination from neutral particles arising from pileup interactions. A muon is accepted if 
the value of the corrected isolation variable is less than 12% of the muon pj. 

A selection based on E™ 1SS provides good separation between signal events and tt background, 
in the case of leptonic decay of the W boson from top decay. The Ey ISS distributions for events 
collected at y/s = 7 and 8 TeV are shown in Fig. [ 2 ] for events with a reconstructed muon pair 
with invariant mass ■ „ > 60 GeV. The background contributions from SM processes are 
















6 


4 Event selection 


Table 2: Event selection: the criteria listed in the upper part of the table are common to the Cl 
and C2 categories, that are then mutually exclusive. 

Common selection 

Single muon trigger pj > 24GeV + isolation + |;/| <2.1 
Event primary vertex |zpv| < 24 cm 

Muon selection 2 opposite-charged muons, 

p T > 24GeV, \t]\ < 2.1, 
track quality cuts, 

\d X y\ < 0.02cm, \d z \ < 0.1cm, 
angular matching with trigger, 
isolation 


£iruss 

E™ 1SS < 35 GeV 

Category Cl 

b tag 

1 or 2 b-tagged jets, 
pif > 20 GeV, |pi et | < 2.4 

Category C2 

No b tag 

Events with no b-tagged jets 


superimposed. For illustration, the expected distribution for signal processes is also shown for 
fflA = 150 GeV and tan /3 = 30. Studies performed using the simulation show that the Ej' ss dis¬ 
tribution for signal events does not vary significantly for different Ota and tan /i assumptions, 
and indicate that the selection E™ lss < 35 GeV provides highest sensitivity for signal at both 
centre-of-mass energies. 

The reconstructed jets are required to have transverse momenta Pj 1 > 20 GeV within the range 
\t] | < 2.4. A multivariate analysis technique is used to remove jets from pileup interactions H5ll . 
Tagging of b quarks in jets relies on the combined secondary-vertex discriminator ||52l , based on 
the reconstruction of the secondary vertex from weakly decaying b hadrons. The discriminant 
frdisc is constructed from tracks and secondary vertex information, and helps to distinguish jets 
containing b, c, or light-flavour hadrons. Jets with an associated frdisc > 0.679 are considered to 
be b tagged. This value represents a good compromise between efficiency to tag b jets in signal 
events from AP («80%) and mistagging probability for light-quark jets (~1%). Figure [3] shows 
the distribution of b L \ lsc in events that satisfy the selection E™ ISS < 35 GeV, for the data collected 
in the two beam energies. For each event, the largest value of b disc is selected. The distribution 
of signal events from the AP process for niA = 150 GeV and tan /l = 30 is superimposed. Jets 
originated from b quark fragmentation tend to be emitted more forward in signal events than 
for tt, thus resulting in a lower observed b-jet multiplicity. For this reason the tt background is 
further suppressed by rejecting events with more than two b-tagged jets, without significantly 
affecting the selection efficiency for signal. 

The events are split into two mutually-exclusive categories. The first category (Cl) contains 
events with at least one jet identified as originated from b-quark fragmentation (b tagged), and 
provides highest sensitivity to AP production channel. Events that do not contain b-tagged 
jets are assigned to category 2 (C2), and provide sensitivity to GF production. The dimuon 
invariant mass distributions for the Cl and C2 categories are shown in Fig. [4] for data and 
simulated events for both centre-of-mass energies. The distributions expected for MSSM Higgs 
bosons with ni\ = 150 GeV and tan /3 = 30, derived from the 1 scenario are also given 
for comparison. A double peak structure around 125 and 150 GeV appears in the C2 category, 
due to the h boson and A+H bosons, respectively. The lower peak is not visible in Cl, as the h 
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Figure 4: The dimuon invariant mass distribution for events that belong to Cl (upper left) and 
C2 category (upper right), for data and simulated events at y/s = 7 TeV. The corresponding 
quantities are shown for y/s = 8 TeV (lower left and lower right). The expected contributions 
to signal assuming the 1 scenario for m A = 150 GeV and tan /l = 30 are displayed for 
comparison. 

production is suppressed in the AP mechanism relative to the GF process. 


5 Signal selection efficiency 

While the calculations for the MSSM cross sections performed in the narrow-width approxima¬ 
tion refer to the on-shell Higgs boson production, at large values of m A and tan /i the convolu¬ 
tion of the larger intrinsic signal widths with the parton distribution functions (PDF) results in 
a non-negligible fraction of signal events produced significantly off-shell. Events with invari¬ 
ant mass significantly smaller than its nominal value have a lower reconstruction efficiency 
than those produced near the mass peak. For consistency, we define signal efficiency as the 
probability for a signal event with the generated invariant mass close to its nominal value to 
be reconstructed and pass all selection requirements of this analysis. The closeness is defined 
using a window of size equal to 3 times the intrinsic signal width (an uncertainty associated 
with this definition is evaluated using a window of 5 times its width, as discussed in Section[7j!. 
With this definition, the product of the MSSM Higgs boson production cross section, luminosity 
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5 Signal selection efficiency 
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Figure 5: Signal efficiency for the AP process at i/s = 8 TeV, shown separately for the three 
(p boson types, (upper left) h, (upper centre) H and (upper right) A, as a function of rn A . The 
corresponding efficiency for the GF production process is shown in the lower row. The con¬ 
tributions from the two event categories Cl and C2 are combined. The results are integrated 
over tan /3, since the efficiency does not strongly depend on this quantity. The band shows the 
change in efficiency due to the limited number of simulated events. 


and signal efficiency provides the normalization for the Higgs boson produced near on-shell. 
The full predicted rate of signal events also contains an additional off-shell contribution, which 
varies with m A and tan /3 and is less than 5% for m A < 250 GeV and tan /5 < 15, and can be as 
large as 15% for rn A = 300 GeV and tan /3 = 30. 

Additional corrections are applied to the signal efficiency to take into account differences be¬ 
tween data and simulation in the muon trigger, reconstruction, and isolation efficiencies. A 
correction is also applied to account for known data-simulation discrepancies in the b tag¬ 
ging efficiency and mistagging probability. The corrections are summarized by a weight factor, 
which is assigned to each signal event. The average of the weight factors computed over all the 
events is very close to one, reflecting the fact that the simulation describes the data with good 
accuracy. 

Figure [ 5 ] shows the signal efficiency at y/s = 8 TeV for AP (top) and GF (bottom) process after 
combining the two event categories Cl and C2. The efficiencies at >/s = 7 TeV are similar. The 
band in the figure represents the variation of the efficiency due to the limited statistics of the 
samples used. The relative amount of AP and GF events in the two event categories varies 
with m a and tan /l, since the production cross sections of the two processes depend on these 
parameters. For example, in the case m A = 150 GeV and tan (1 = 30, more than 90% of the 
signal events in Cl would be from AP production, and about 60% in C2. For m A = 150 GeV 
and tan (5 = 5, where the GF contribution becomes more relevant, the content of AP events 
would be 60% in Cl and only 15% in C2. 
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6 Fit procedure 

The procedure described below is applied separately to Cl and C2 events. The event selection 
criteria are applied to the simulated samples listed in Table [l] For each sample, and for each of 
the three cp bosons, the invariant mass distribution of the events that pass the event selection 
is approximated with a Breit-Wigner function convolved with a Gaussian, that accounts for 
detector resolution. This analytical expression provides a good description of the signal shape 
for all the m A and tan f> values. The three functions are denoted F^, Fh, and Fa, and contain 
the mass and width of the Breit-Wigner and the width of the Gaussian as free parameters. The 
function F S i g represents the expected signal yield, and it is a linear combination of the three 
functions described above: 

F sig = w h F h + w H F h + w A F A , (2) 

where w h , ?o H , and w,\, are the number of events containing h, H, and A bosons, respectively, 
calculated according to their expected production cross sections. An example of this procedure 
is shown in Fig. [6](left) for m A = 150 GeV and tan f> = 30. The highest peak represents the 
superposition of the contributions from H and A bosons, that in this case are almost degenerate 
in mass. 

CMS Simulation 19.3 fb 1 (8 TeV) CMS 19.3 fb 1 (8 TeV) 




Figure 6: Invariant mass distribution of the expected signal for m A = 150 GeV and tan f = 30 
(left), and an example of the fit to the data at i/s = 8 TeV including the same signal assumption 
(right). The distribution represents the expected number of events for an integrated luminosity 
of 19.3 fb -1 . For each plot the pull of the fit as a function of the dimuon invariant mass is shown. 


Since the Drell-Yan muon pair production is the dominant background process, it is modeled 
by a Breit-Wigner function plus a photon-exchange term, which is proportional to l/m^+ . 

Defining m = mp+p-, the function F bkg becomes: 



fz _1_ 

N r m (m - m z ) 2 + r f 


(Hz)i 

N n°rm ? „2 ' 


0 ) 


where e Am describes the effects of the PDF, and the /V ( norm terms correspond to the integral of the 
corresponding functions in the chosen mass range. The quantity fz represents the contribution 
of the Breit-Wigner term relative to the photon-exchange term. The quantities A and fz are 
free parameters of the fit. The parameters Vz and niz are determined separately for the Cl and 
the C2 events from a fit to the in,, f , distribution in the mass range of the Z boson between 80 
and 120 GeV. The fit provides the effective values of such quantities, that include detector and 
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7 Systematic uncertainties 


resolution effects for each set of data. Their values are used in F^g and are kept constant in the 
fit. 

A linear combination of the two functions for the expected signal and background is then used 
in an unbinned likelihood fit to the data: 


Tfit = (1 ~ /bkg) fsig + /bkg fbkg- (4) 

The parameters that describe the signal are determined in the fit of the simulated signal to 
Eq. §, for each pair of m\ and tan / values. Subsequently, they are fixed in Ff it , where the 
free parameters are the quantities A, fz, and /bkg- The fraction of signal events is defined as 
/sig = (1 — /bkg)- The data are fitted to Ffi t in the mass range from 115 to 300 GeV for each point 
in the m A and tan / parameter space. As an example, the fit to the data of C2 at = 8 TeV is 
illustrated in Fig. [6](right), assuming a signal with m A = 150 GeV and tan / — 30. 


7 Systematic uncertainties 

The following sources of systematic uncertainties are taken into account, and the impact of one 
standard deviation change is reported in terms of a variation in the nominal signal efficiency 
defined in Section |5] 

The limited number of simulated events introduces an uncertainty in the signal selection ef¬ 
ficiency that is at most 2.0%. The muon trigger, reconstruction, identification, and isolation 
efficiencies are determined from data using a tag-and-probe technique lf38| . The uncertainty in 
the trigger efficiency correction is 0.5%, whereas 1.0% is assigned to the combination of uncer¬ 
tainties in muon reconstruction and identification, as well as on isolation efficiencies. 

A systematic uncertainty in the pileup multiplicity is evaluated by changing the total cross 
section for inelastic pp collisions in simulation. The corresponding uncertainty on the signal 
efficiency is at most 0.8% in both categories. 

The event fractions in the two categories depend on the b tagging efficiency and the mistag- 
ging probability. The uncertainty in the b tagging efficiency is estimated by comparing data 
and simulated events with samples of enriched b quark content and different topologies, as 
described in Ref. If52l . The uncertainty in the efficiency to detect b jets is about 3.0%. Similarly, 
the uncertainty in the mistagging rate is about 10%. Their overall contribution to the selection 
efficiency is weighted by the fraction of AP and GF events that are expected in each event cate¬ 
gory, which depends on »za and tan /. The largest overall uncertainty is 3.0% for Cl, and 0.4% 
for C2 events. 

The jet energy scale uncertainty is estimated by smearing the jet momentum by a factor depend¬ 
ing on pj and ;/ of each jet, as described in Ref. |37j. The effect on signal selection efficiency 
is 4.0% for events that belong to the Cl and 0.5% for the C2 categories, at \fs = 8 TeV. For 
y/s = 7 TeV the corresponding numbers are 3.8% and 0.6%. The uncertainty in the £/" ss scale 
and resolution is estimated through comparisons between data and simulation 11531 541- The 
effect on the signal selection efficiency is 3.0% and 2.0%, the same for both categories, for the 
sample with yjs = 8 and 7 TeV, respectively. The uncertainty in the integrated luminosity is 
2.6% and 2.2% at i/s = 8 and 7TeV, respectively Il55ll56). 

Uncertainties due to the choice of PDF set affect the signal efficiency, and are studied using the 
PDF4LHC It57l prescription. The renormalization and factorization scales in the calculations 
and their changes are summarized in Refs. lfl6l - [T8l . The effect on the signal selection efficiency 
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varies from 1.0% to 3.0% over the m ,\ and tan fi parameter space. The choice of 3.0% is taken as 
the systematic uncertainty 

The efficiency is determined for events with generated mass values within a window of a factor 
of 3 of the intrinsic width of the Higgs boson, as described in Section [5j The difference relative 
to the efficiency obtained using a cutoff of a factor of 5 of the intrinsic width is assigned as a 
systematic uncertainty The uncertainty is between 1% to 3% for the Cl and 1% to 5% for the 
C2 categories. 

Table [3] lists the systematic uncertainties that affect the determination of signal efficiency. The 
impact of these systematic uncertainties on the exclusion limits that will be presented in Sec¬ 
tion [8|is negligible compared to the statistical uncertainty. All the systematic uncertainties in 
Tableppre correlated for the v^s = 7TeV and 8 TeV data, with the exception of the uncertainties 
related to the limited MC statistics and the integrated luminosity 

The uncertainties in the MSSM cross sections depend on m A , tan /3, and the scenario, and are 
provided by the LHC Higgs Cross Section Working Group Iil6l - (l8l . The signal events are gen¬ 
erated using pythia, assuming the parameters of the m“ ax scenario, as discussed in Section [3] 
The different benchmarks are expected to affect the production cross section, but not the kine¬ 
matic properties of the events related to Higgs boson production and decay. To check this 
assumption, events are generated with PYTHIA using the parameters for the m ™ oA+ , m™ od ~, 
light stop and light stau benchmarks, assuming rn A = 150 GeV and tan /3 = 20. The events are 
generated for both the GF and the AP mechanisms, and the Higgs boson pj and the E™ 1SS of 
the events are compared at generator level for the various benchmark scenarios. No significant 
differences are observed in the distributions of these quantities. 

Table 3: Sources of systematic uncertainties for Cl and C2 event categories that affect the signal 
efficiency at y/s = 8 TeV. They are expressed in terms of relative signal selection efficiency 
When the systematic uncertainty at i/s = 7 TeV differs from y/s = 8 TeV, the corresponding 
value is quoted in parenthesis. 


Source 

Systematic uncertainty (%) 

Cl 

C2 

MC statistics 

2.0 

2.0 

Trigger efficiency 

0.5 

0.5 

Muon efficiency 

1.0 

1.0 

Muon isolation 

1.0 

1.0 

Pileup 

0.8 

0.8 

b tagging 

3.0 

0.4 

Jet energy scale 

4.0 (3.8) 

0.5 (0.6) 

£miss 

3.0 (2.0) 

3.0 (2.0) 

Integrated luminosity 

2.6 (2.2) 

2.6 (2.2) 

PDFs 

3.0 

3.0 

Width correction 

1-3 

1-5 


Since the number of background events is determined through a fit to the data, an additional 
systematic uncertainty arises from the possibility that the background parametrization may 
not adequately describe the data as a function of the dimuon invariant mass. A method sim¬ 
ilar to that described in Ref. IflO] is used to evaluate the effect, by estimating the uncertainty 
through the bias in terms of the number of signal events that are found when fitting the signal 
+ background model (as described in Section [6j to pseudo-data generated for different alter¬ 
native background models. Such alternative background parametrizations include Bernstein 
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8 Results 


polynomials and combinations of Voigtian and exponential functions. Bias estimates are per¬ 
formed for mass points between m a = 115 and 300 GeV. For each m,\ value, the largest bias 
among the tested functions is taken as the resulting uncertainty. The bias is implemented as a 
floating additive contribution to the number of signal events, constrained by a Gaussian prob¬ 
ability density with mean of zero and width set to the systematic uncertainty. The width of the 
Gaussian is the largest systematic uncertainty, and the effect is to increase the expected limit on 
the presence of a signal by 20% in the region near via = 120 GeV and by about 10% at larger 
mass values. 

In the mass range between 115 and 300 GeV, that is relevant for this analysis, the mass reso¬ 
lution is estimated to be between 1.2 and 4 GeV. Uncertainties in the muon momentum deter¬ 
mination can affect the invariant mass measurement, and have been carefully studied in data 
and simulation [[38|. The dimuon invariant mass resolution for masses above the Z peak has 
been previously studied in the search for a SM Higgs decaying to a dimuon pair 11581 - The mass 
resolution determined from data at the Z mass value is 1 GeV, in excellent agreement with the 
prediction from simulation. This value is consistent with the mass resolution of 1.2 GeV that 
we estimate from simulation for a mass of 115 GeV, that corresponds to the lower edge of the 
Higgs mass range considered in this analysis. 

The overall capability of the analysis to detect the presence of a signal is verified by introduc¬ 
ing a hypothetical simulated signal in the data using the shape parametrization discussed in 
Section [6] The average measured number of signal events is found to be within 1.3% of the 
injected signal for the Cl category, and within 4.3% for the C2 category. These differences are 
assigned as systematic uncertainties. 


8 Results 

No evidence of MSSM Higgs bosons production is observed in the mass range between 115 
and 300 GeV, where the analysis has been performed. Upper limits at 95% confidence level 
(CL) on the parameter tan j6 are computed using the CL S method [59, 601, which is a modified 
frequentist criterion, and are presented as a function of m/\. Systematic uncertainties are incor¬ 
porated as nuisance parameters and treated according to the frequentist paradigm |6l1 . The 
results are obtained from a combination of both event categories and centre-of-mass energies. 
For each value of ui f \, the value of tan /3 at which the CL exceeds 95% is chosen to define the 
exclusion limit on that m ,\. This is performed for all the rn,\ values and the results are shown in 
Fig. 0 These results are obtained within the m™ od+ scenario. The observed upper limits range 
from tan ft of about 15 in the Iow-twa region, to above 40 at m ; \ = 300 GeV. For larger values of 
niA the uncertainty on the tan /3 upper limit becomes large, exceeding tan /3 = 50, for which the 
MSSM cross-section predictions are not reliable. 

A comparison with the results obtained for the m ™ d ~, m™ ax , light stop and light stau scenarios 
is also performed. The exclusion limits computed within these other benchmark models are all 
very similar. For any value of niA, the quantity A tan j6 = tan /3 (); m 0 d — tan /^scenario represents the 
difference of the tan (5 values at which the 95% CL limit is determined if an alternative scenario 
is used. Figure [8] shows the quantity A tan p as a function of rn,\ for all the tested scenarios. For 
most niA values, the 95% CL limits on tan f ] > computed within a given scenario differ by less 
than one unit from the results obtained within the m™ od+ scenario. 

Limits on the production cross section times decay branching fraction a B(cp —t }i~ pi ) for a 
generic single neutral boson (p are determined. In this model independent analysis no assump¬ 
tion is made on the cross section, mass, and width of the (p bosons, which is sought as a single 
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Figure 7: The 95% CL upper limit on tan /3 as a function of rn,\, after combining the data from the 
two event categories at the two centre-of-mass energies (7 and 8TeV). The results are obtained 
in the framework of the m™ od+ benchmark scenario. 

resonance with mass m^. The analysis is performed assuming the narrow width approxima¬ 
tion, for which the intrinsic width of the signal is smaller than the invariant mass resolution. 
For this purpose the simulated signal of the A boson for the case tan ^ = 10 is used as a template 
to compute the detection efficiency for a generic (p boson decaying to a muon pair. The single 
(p boson is assumed to be produced entirely either via the AP or the GF process, and the search 
for a single resonance with mass is performed. The 95% CL exclusion on aB((p —? y. + }C) 
is determined as a function of m^, separately for the two production mechanisms. The combi¬ 
nation of events belonging to Cl and C2 is shown in Fig. [9j assuming the (p boson is produced 
either via the AP or the GF process. Only data collected at i/s = 8 TeV are used, as they pro¬ 
vide a better sensitivity because of the higher luminosity. In addition, since the (p production 
cross section depends on the centre-of-mass energy, a combination with the 7 TeV results would 
introduce a model dependence in the description of the cross section evolution with energy. 


9 Summary 

A search has been performed for neutral MSSM Higgs bosons decaying to ji 1 }t from pp col¬ 
lisions collected with the CMS experiment at \Js = 7 and 8 TeV, corresponding to integrated 
luminosities of 5.1 and 19.3 ftC 1 , respectively. The analysis is sensitive to Higgs boson produc¬ 
tion via gluon fusion, and via association with a bb quark pair. The results of the search, which 
has been performed in the mass range between 115 and 300 GeV, are presented in the m™ od+ 
framework of the MSSM. With no evidence for MSSM Higgs boson production, this analysis 
excludes at 95% CL values of tan /3 larger than 40 for Higgs boson masses up to 300 GeV. Com¬ 
parisons with m™ od m™ ax , light stop, and light stau scenarios are also presented, and offer 
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9 Summary 



Figure 8: Comparison of the 95% CL exclusion limits on tan ft obtained within MSSM bench¬ 
mark models, as a function of m\. The quantity A tan ft = tan ft m mod+ — tan ft scena no represents 
the difference in tan ft at which the 95% CL limit is obtained for alternative scenarios. 




Figure 9: The 95% CL limit on the product of the cross section and the decay branching fraction 
to two muons as a function of niq, obtained from a model independent analysis of the data. The 
results refer to (left) b quark associated and (right) gluon fusion production, obtained using 
data collected at y/s = 8 TeV. 


very similar results relative to the m™ od+ benchmark. Limits are determined on the product 
of the cross section and branching fraction crB(cp —> }t + ^ ~) for a generic neutral boson cp at 
y/s = 8 TeV, without any assumptions on the MSSM parameters. In this case the (p boson is 
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assumed to be produced either in association with a bb quark pair or directly through gluon 
fusion, and sought as a single resonance with mass m ( p. Exclusion limits are in the mass region 
from 115 to 500GeV. For nip = 500GeV, values aB{(p —> y + }r) > 4fb are excluded at 95% CL 
for both production mechanisms. These are the most stringent results in the dimuon channel 
to date. 
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